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Nonperturbative, in inverse Thouless conductance , corrections to distributions of level veloc- 
ities and level curvatures in quasi-one-dimensional disordered conductors with a topology of a ring 
subject to a constant vector potential are studied within the framework of the instanton approxima- 
tion of nonlinear cr-model. It is demonstrated that a global character of the perturbation reveals the 
universal features of the level dynamics. The universality shows up in the form of weak topological 
oscillations of the magnitude ~ covering the main bodies of the densities of level velocities and 
level curvatures. The period of discovered universal oscillations does not depend on microscopic pa- 
rameters of conductor, and is only determined by the global symmetries of the Hamiltonian before 
and after the perturbation was applied. We predict the period of topological oscillations to be 4/-jt^ 
for the distribution function of level curvatures in orthogonal symmetry class, and VS/n for the 
distribution of level velocities in unitary and symplectic symmetry classes. 
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I. INTRODUCTION 



Parametric level statistics describes the response of the 
spectrum of complex chaotic systems to an exter- 

nal perturbation, being a measure of the sensitivity of 
the individual energy level to a change of some exter- 
nal parameter x. The part of the parameter x can be 
played by either the external electric or magnetic field, 
or the 'strength' of a background potential. It may be 
any other relevant parameter of the Hamiltonian. 

In the series of papers it was shown that a sys- 

tem whose spectrum follows closely the universal fluctu- 
ations predicted by the Wigner-Dyson random matrix 
theory Q (RMT) should also exhibit a universal para- 
metric behavior. This conclusion was reached by ana- 
lyzing the correlations of the single electron level den- 
sities in a disordered metallic sample with the topology 
of a ring pierced by an Aharonov-Bohm magnetic flux 
if. In this particular problem, the parametric correla- 
tions take a universal form involving the rescaled pa- 
rameter oc gip'^, with the dimensionless conductance 
g = Ec//^, the ratio of the Thouless energy and the 
mean level spacing. Numerical simulations |2j have sup- 
ported the point that the universal character of para- 
metric level statistics extends to a wider class of chaotic 
systems without disorder (chaotic billiards) whose Hamil- 
tonian depends on some external parameter x. In such 
systems, the spectral fluctuations taken at different val- 
ues of X become system-independent after the rescaling, 
x t-^ X, which involves solely the 'generalized' dimen- 
sionless conductance g cx (^[dEn{x) / dx]^^Q) . [Here, 
the angular brackets (...) stand for disorder averaging.] 
Further investigations have established the univer- 
sality of other parametric statistics such as the distribu- 



tion functions P^{V) = A{J2nS{V - Vn)S{E - £;„(0))) 
andPciK) ^ A{J^^6{K~Kn)S{E^En{Q))) of the level 
velocities Vn = A~-^[dEn{x) / dx]x^Q and the level curva- 
tures Kn = A~^[d'^En{x)/dx^]x^o- These two measures 
of level dynamics turn out to be universal after the rescal- 
ing V = V{V^y/^ and K = K{\Kn\). Both distributions 



Pv{V) cx 



/3=1, 



exp{-^V2}, /3-2,4, 



and 



PciK) cx (1 + 
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(1.1a) 



(1.1b) 



appear to depend only on the Dyson index /3 = 1,2 
and 4. In turn, the averages (V^) and (l-ftTnl) are deeply 
connected, via the Edwards-Thouless relation |9| , to the 
average system conductance if the parameter x repre- 
sents the Aharonov-Bohm flux. The latter circumstance 
makes it possible to use these rather abstract measures of 
level dynamics for distinguishing between the extended 
and the localized phases in disordered systems. 

The above universality (that emerges after appropriate 
rescaling) implies that the resulting parametric correla- 
tions do not depend on details of the system and the 
perturbation. Instead, they are only determined by the 
fundamental (orthogonal, unitary or symplectic) symme- 
tries of the system before and after it was perturbed. For 
disordered conductors, this universal regime is observed 
in the limit of the inflnite dimensionless conductance, 
g = £'c/A — !■ oo. In quantum chaotic systems, the same 
criterion pd| applies provided one identifies the Thouless 
energy with the first nonzero mode 71 in the spectrum 
of the Perron-Frobenius operator of the corresponding 
classical system. It should be stressed that the extended, 
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structureless electron eigenf unctions which do not corre- 
late with the fluctuations of the electron energy levels 
are the underlying physical reasons of the universality 
phenomenon. 

One-dimensional conductors in the strongly localized 
regime constitute another extreme situation, in which the 
universality of parametric statistics is broken completely. 
Recent studies have shown that, in the regime of 

strong localization, the distribution function of level cur- 
vatures has nothing in common with the universal law 
Eq. (1.1b). Instead, it roughly follows the log-normal 



law which explicitly contains such system-dependent pa- 
rameters as the circumference L of a one-dimensional ring 
and the electron mean free path I (which is of order of 
the localization length ^ for a one-channel conductor). 

Disordered conductors with large but finite values of 
g, which are in between of the two extreme situations 
described above, offer a fertile field for tracing the mech- 
anism of breaking the parametric universality. While, in 
this case, the electron wave functions are still extended, 
they do have a pronounced internal structure (showing 
up in their long-range correlations [^,^) that cannot 
already be neglected. Finite-g corrections to paramet- 
ric level statistics beyond the Wigner-Dyson RMT have 
been a focus of a number of studies |16-18|. All of them 
have emphasized a nonuniversal character of the level 
dynamics manifesting itself already in the form of the 
perturbative corrections ||l^] to the level curvature dis- 
tribution which becomes sensitive to both the spatial di- 
mensionality d of the disordered conductor and the origin 
of the external perturbation. The latter determines the 
sign of the correction which depends on the presence or 
absence of the global gauge invariance (that is, on the 
topological origin of the perturbation) and is different 
for the T-breaking perturbations in the form of random 
magnetic fluxes that act locally and for the T-breaking 
perturbations represented by a constant vector potential 
(that is equivalent to a global phase twist in the bound- 
ary conditions). Nonanalyticity |^ of the level curvature 
distribution in the latter case in the low-dimensional con- 
ductors {d— 1,2) was also addressed. 

In the present paper we further examine the level dy- 
namics in quasi-one-dimensional disordered conductors 
in the presence of the Aharonov-Bohm magnetic flux to 
demonstrate that, for finite conductance g, a new kind of 
topological universality emerges, for which the existence 
of the global gauge invariance is crucial. This is a non- 
perturbative in g^^ effect. To be precise, we argue that 
the ring topology and the global character of the pertur- 
bation are responsible for appearance of the weak oscil- 
lations covering the main bodies, Eqs. (1.1a) and (l.lt), 
of the distributions of the level velocities and the level 
curvatures. While the magnitude of the oscillations is 
system dependent (being of order of ^ e"^), their period 
appears to be independent of microscopic parameters of 
the system and the perturbation: It is entirely deter- 



mined by the global symmetries - orthogonal, unitary or 
symplectic - of the system before and after the pertur- 
bation was applied, and is different for level velocities 
and level curvatures. We also stress that, contrary to the 



Wigner-Dyson universality [Eqs. (1.1a) and (1.1b)], that 
arises after rescaling V i-^ V and K ^-^ K , no rescaling 
is needed to establish the universality of topological os- 
cillations: Their period is universal and parameter inde- 
pendent in genuine variables V and K . In the absence of 
the global gauge invariance, the above universality does 
not show up at all. 

To appreciate the difference between the local and the 
global perturbations on the formal level, it is instructive 
to appeal to the random matrix formulation of the prob- 
lem, replacing the microscopic Hamiltonian by banded 
random matrices which are known |po| to describe (in a 
certain thermodynamic limit) the physics of disordered 
quasi-one-dimensional conductors. For definiteness, the 
orthogonal symmetry is supposed to be respected in the 
unperturbed conductor, while the perturbation is as- 
sumed to drive the system symmetry toward the unitary 
one. 

(i) For the local perturbation of the strength x, the 
perturbed quasi-one-dimensional disordered system with 
N sites can be modeled by the Nx N parametric random 
matrix of the Pandey-Mehta type |2^] 



i/'°'=^'(a;) = H' 



ixH' 



(1.2) 



where and H"^ are N x N real symmetric and real 
antisymmetric, statistically independent banded random 
matrices, respectively. Notice that such a decomposition 
imposes no constraints for the real eigenvalues of the ma- 
trix H^°'^^^{x) as functions of the variable x. 

(ii) This contrasts the case of the global perturbation 
applied to quasi-onc-dimensional conductors with a ring 
topology which should explicitly be incorporated |^ to 
both unperturbed and perturbed matrix Hamiltonians. 
An unperturbed (real symmetric) matrix H{0) must ac- 
count for a periodic electron motion along a discrete ring 
of N sites (see Fig. 1). 




FIG. 1. A discrete ring of N sites pierced by a (dimension- 
less) Aharonov-Bohm flux if. Connected geometry of a ring 
permits the electron hopping from the A'^-th site to b adjacent 
sites to the left and b sites to the right. 
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FIG. 2. Schematic representation of the structure of the 
periodic banded matrix that contains three (shaded) regions 
with nonvanishing entries: The band of the width 2b along 
the main diagonal, the upper right and the lower left corners 
(both of the width b). In the presence of the Aharonov-Bohm 
flux, the c orne r entries are multiplied by the phase factors 



Eq. (0. 



If the particle located at some site may reach only 
b neighbors to the right and b neighbors to the left 
(6 < [N/2]), it is easy to see that the unperturbed 'peri- 
odic' N X N matrix H{0) will contain three regions with 
nonvanishing entries: a band \^J. — i^l < b along the main 
diagonal, and the upper right and the lower left corners, 
\^j.-iy\> N -b (see Fig. 2). 

Let us now apply a small static magnetic flux $ = 
((^/27r)$0j 'i'o — hc/e being the flux quantum. Each 
time the electron hops from the site fi to the site /i + 1 , it 
ascribes the phase (f/N. As a result, the transition am- 
plitude between the sites fi and i/ gets modifled by the 
phase factor e***"", Hni, Hn,ye^^^'", with 



sgn{fi -iy)[l-\i^- fi\/N]ip if\fi-iy\> N ~b. 

(1.3) 



Passing on to the basis, in which the eigenvectors of the 
perturbed matrix are transformed as tpk '-^ VfeC*'''^^'''^, 
we arrive at the perturbed periodic matrix of the form 

r H^,{0) if\fi-iy\<b, 
[ iJ^^(0)e*'^ if^i-u>N-b. 

The perturbed periodic matrix iJgi°bai obtained in this 
way possesses a manifestly different structure as com- 
pared to the structure of H^°'^^^^ hereby reflecting the 
presence of the global gauge invariance in the former case. 

The structure of Eq. (L4) dictates the 27r-periodicity 
of the eigenvalues En{ip) of H^'^°^'^'^{ip), 



(1.5) 



Equation (^^) suggests the Fourier analysis to be an ap- 
propriate language for description of level dynamics in 
the case of the global perturbation: 



A cos(to^). 



(1.6) 



As a matter of fact, Eq. (1.6) allows one to re- 
late the problem of evaluating the Fourier transform 
Pc(A) — J dK exp{—iKX}Pc{K) of the curvature dis- 
tribution function to an effective statistical mechanical 
problem. To see this, we notice that the level curvature 
Kn can be expressed in terms of the amplitudes 



as 



1 d^En 



A 9(^2 



E 

rn— 



(1.7) 



so that the Fourier transform Pc(A) of the curvature dis- 
tribution function is 

P,(A) = A exp I ^X m\^^^ \ 5{E ~ E^m) • 

\ n I m=0 J / 

(1.8) 

Introducing the joint probability distribution function for 
the amplitudes in the form 

P,({a^}) = A 5{E - £;„(0)) \{ 5(a™ - a^")) \ , 

\ n m— / 

(1.9) 

we conclude that the free energy Fc{X) = — InPc(A) of 
the effective statistical mechanical model equals 



M 



F,{\) = - lim In / TT da™Pg({a™}) 

X exp < i\ ^ m^a„i > . (1-10) 



k m=0 ) 

Here, Pg({am}) describes the interaction of fictitious 
'particles' subject to external 'field' i\. [Similar repre- 
sentation can be obtained for the distribution function of 
the level velocities Vn for /3 = 2,4. In this case, cos(m(/3) 
is to be replaced by sin(m(y9) in the Fourier expansion 
Eq. (1.6).] It is important to emphasize that no such 
mapping to an effective statistical mechanical problem is 
available for the case of the local perturbation, Eq. (|l.2|). 

In the rest of the paper, we argue, by using the instan- 
ton approximation [psljiol of nonlinear cr-model, that in 
the li mit g ^ 1 the effective statistical mechanical model 
Eq. ( 1.10 ) exhibits a second-order phase transition as 
the external field i\ varies, with the free energy Fc{\) 
displaying a discontinuity of its second derivative at the 
'critical' point determined by the value of the critical 
parameter Ac = 7r'^/2. It turns out, that the value of 
the critical field iX^ is universal: It does not depend on 
microscopic parameters of the problem but is entirely 
determined by the global character of the perturbation, 
Eq. ( |1.4D , that breaks the orthogonal symmetry toward 
the unitary one. This finding, being translated to the 
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space of genuine curvature variable K, implies appear- 
ance of the universal oscillations with the universal pe- 
riod 6K — 27r/Ac — 4/7r^ and the system-dependent am- 
plitude ~ on the main body Eq. ( |l.lb| ) of the level 
curvature distribution. Similar analysis of the distribu- 
tion functions of the level velocities {at P — 2 and 4) 
reveals the existence of the second-order phase transition 
in associated effective statistical mechanical problem as 
well. In the latter problem, the value of the critical exter- 
nal field iXc is also found to be universal though different: 
Ac = 27r^/-\/3. As a result, the main body Eq. (1.1a) 
of the distribution function of the level velocities gets 
dressed by the weak oscillations of the same magnitude 
(~ e~^) with the universal period 6V = 27r/Ac = -s/S/tt. 
We stress that the above universality has a clear topolog- 
ical origin: A global character of the perturbation (con- 
stant vector potential) and a ring geometry of a quasi- 
one-dimensional conductor are necessary conditions for 
its existence. 

The paper is organized as follows. Section II is devoted 
to a nonperturbative analysis of the level curvature dis- 
tribution function associated with the global T-breaking 
perturbation over the orthogonal ensemble. The distri- 
bution of the level velocities in ensemble with the uni- 
tary symmetry is analyzed in Section III. In Section IV 
a relationship between the ring geometry of a conductor, 
the boundary conditions in associated nonlinear u-model 
(in its quasi-one-dimensional and ergodic limits) and the 
discovered topological oscillations is discussed. Section 
V contains conclusions. The details concerning the dis- 
tribution function of the level velocities in the case of 
symplectic symmetry are collected in Appendix. 



II. DISTRIBUTION OF LEVEL CURVATURES 

AT /3 = 1 



Let us consider the electron of mass m moving in a 
disordered conductor with a topology of a quasi-one- 
dimensional ring of the circumference L. In the pres- 
ence of the magnetic flux $ that pierces the ring, the 
one-particle Hamiltonian takes the form 



H 



2;^K + z^^ 



Vir), 



(2.1) 



where ip — 27r$/<I>o, 'i'o = hc/e is the flux quantum, 
and eg is the unit vector along the azimuthal direction. 
The potential ^(r) describes the effect of disorder and 
is supposed to be the 5-correlated random Gaussian pro- 
cess with {V{r)V{r')) = {n/2'KVT)5{Y - r'). Here, r is 
the mean time between collisions and v is the density 
of states. The (dimensionless) magnetic flux Lp in the 
Hamiltonian breaks the time-reversal invariance which 
has been respected in the absence of the perturbation. 



In what follows we are interested in the distribution 
function Pc (if ) = /:^{Y.n^iK - Kn)5{E - En{Q))) of the 
level curvatures 



2 E^{^) - 
Kn^ -r lim 5- 



(2.2) 



fluctuating, for a given level, over the ensemble of disor- 
der realizations. These fluctuations are due to fluctua- 
tions of both the single electron eigenfunctions and the 
electron eigenlevels in an unperturbed system. [Here, 
A is the mean level spacing.] Distribution function of 
level curvatures can be expressed through the paramet- 
ric two-level correlation function R{uj,(p) — A'^{i^{E + 
uj, tf)v{E, 0)) by means of the limiting procedure p7[| 



Pc{K) 



lim 



R uj = 



(2.3) 



similar to the one suggested in Ref. for the distribu- 
tion of the level velocities. The relation Eq. ( pT^ ) is based 
on the fact that the level velocity vanishes identically in 
ensembles with the weakly broken orthogonal symmetry, 
Eq. ( pTTal ). 

The parametric two-level correlation function i?(w, i^), 
and hence the distribution function Pc{K) of the level 
curvatures, can be represented in the form of the func- 
tional integral over the 8x8 supermatrix field Q by using 
the Efetov's nonlinear cr-model [|To| . The latter statistics 
is given by [pljri 



Pc{K) ^ lim (^2Re f VQp4Q] exp{-F[Q; K, 



(2.4) 



Here, Pc[Q] is some preexponent irrelevant for our further 
consideration that holds with the exponential accuracy. 
The free energy functional F[Q; K, ip] reads: 



F[Q; K,^]^^J dxStr {O^Q - [f, Q]f 



dxStr (AQ) . 



(2.5) 



In the above formulas, x = r/L G (—1/2, +1/2) is 
the dimensionless coordinate measured along the ring, 
g = TiD/L^A is the dimensionless conductance, g 3> 1. 

For the ring topology the supermatrix field Q{x) must 
obey the periodic boundary conditions. Its symme- 
try properties are specified in Ref. |Q, so that Q{x) 
spans the three subspaces: retarded-advanced (R-A), 
boson-fermion (B-F), and the subspace related to the 
time reversal (T-R). The supertrace, Str, is defined as 
StrM = TiMp F - TtMb b ■ Other matrices entering Eq. 
(PI) are 



Is + A, 



-T3, 



(2.6) 



with 



4 



and 



fs = diag {az,crz,(j^,az)^_^ , 
A = diag (I2, 12, ~l2, -12)r_a ' 



1 



-1 



(2.7) 



(2.8) 



being the z-Pauh matrix in the T-R space. The matrix 
A describes breaking the symmetry between the retarded 
and the advanced subspaces, while the matrix breaks 
the time reversal symmetry. 

Notice that the operator dx — i(p[T,...] in Eq. (2.5) 
implies the presence of the global gauge invariance due 
to the global character of the applied T-breaking pertur- 
bation represented by the Aharonov-Bohm flux (p. This 
operator gets destroyed if the perturbation has a local 
origin, that is the case of random magnetic fluxes or mag- 
netic impurities. In this local case, the term linear in (p 
is absent in Eq. (2.5). 

It was realized in Refs. that the level curvature K 
enters the free energy functional Eq. (2.5) in the same 



way as does uj in the problem of the long-time current 
relaxation This observation suggests that the in- 

stanton approximation of nonlinear cr-model should be 
applicable to the evaluation of the Fourier transform, Eq. 



(1.8), of the level curvature distribution: 



Pc{X) = lim ifi'^Re 

X [ dK [ VQp,[Q]e 



-{F[Q:K,ip]+2K\} 

(2.9) 



Equation (2.9) holds with the exponential accuracy, and 
is valid whenever the functional -F[(5ins; ^sp, 0] +iKspX is 
large enough. [Hereafter, we consider only positive values 
of A because of the evenness of Pc{K).] 

An important remark is appropriate here. In accor- 
dance with the results of Ref. j23| one can expect that at 
large enough A the space-independent configurations of 
Q{x) are energetically unfavorable: There are essentially 
space-dependent, periodic superfields Qinsi^) which min- 
imize the functional in question. In contrast, at relatively 
small A we expect that the minimum of F[Q; K, ip] + iKX 
in the limit (/3 — s- is reached on spatially independent 
configurations of Q{x). This formal argumentation sug- 
gests the existence of some 'critical' point A = Ac at 
which a transition between the two solutions occurs. Be- 
low, this scenario will be confirmed by explicit calcula- 
tions. It turns out that at the critical point the func- 
tional F[Qins', Ksp, 0] + iKspX and its first derivative with 
respect to A are continuous functions of A, and only the 
second derivative exhibits a discontinuity with immedi- 
ate implications for the distribution function Pc{K) itself 
to be discussed later on. 



Two observations are to be made before we proceed 
with the computation of Pc(A). First, we note that due 
to the limit 1^ ^ to be taken in Eq. ( |2.9| ) the only 
nonzero contribution to Pc{X) can arise from the entries 
in the noncompact boson-boson sector of the superma- 
trix Q{x) because they are allowed to take arbitrary 
large values. Second, the integration over the Grass- 
mann entries oiQ(x) modifies only the preexponent Pc[Q] 
in Eq. ( |2.9D and, therefore, they do not contribute to 
Pc(A) = exp{— Fc(A)} if the latter is evaluated with the 
exponential accuracy, provided Fc{X) :s> 1. For these two 
reasons, we may retain only the commuting entries in the 
noncompact sector of Q{x) in Eq. (^!g|). Hence, the su- 
permatrix Q may be thought of as a conventional 4x4 
matrix Qb{x). 

Using the Efetov's parameterization [|o| for the or- 
thogonal symmetry class, we obtain 



Qb{x) = U{x) . . ^ X ] U {x), 



where 



and 



U{x) 



e^x'"- 



(2.10) 
(2.11) 



R-A 



cosOb = cosh^i cosh02l2 + sinh(?i sinh(?20'j;, (2.12a) 
siwOb = i sinh^i cosh 6*212 + * cosh ^1 sinh^^20■a;• (2.12b) 
Here, 



1 

1 



(2.13) 



is the x-Pauli matri x in the T-R space; z-Pauli matrix 
was specified in Eq. (2 



In the equations above, 0, x G (OjStt) and 6*1,^2 G 
(0, +00) are spatially dependent, periodic functions of 
the variable x € (—1/2,-1-1/2). Obviously, this periodic- 
ity reflects the ring topology of the disordered conductor 
in question. Substituting Qb{x) in the parameteriza- 
( 2.12 ) into Eq. (2_^), and varying 
'iKX with respect to the func- 



tion Eqs. ( p.lOD 
the functional F[Q\K,(p\^ 
tions ^1(2;), 92{x), 4i{x) and xi^) find ^ set of four 
saddle-point equations which are satisfied by the choice 
9i{x) — 92{x) — 6{x) and 4>{x) = const. The functions 
9{x) and xi^) obey the following saddle-point equations: 



-^[26»(a;)] +sinh[26i(a;) 



—Xix) - ip 
dx 



d 
dx 



^X(a;) - ip ] sinh^ 9{x) 



= 0, 



= 0, 



(2.14) 



(2.15) 
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where we have introduced k = iKsp/2g. Also, variation 
over K yields the first self-consistency equation 



dx cosh^ 6{x) — 2A. 



(2.16) 



Fc{X) = lim lugidx 



sinh 



Notice that Eq. (2.17) is justified as long as Fc{\) ^ 1. 

In order to perform the limit iy9 — s- 0, we introduce two 
new functions v{x) — x{x)/(p and 9{x) = 9{x) + Ini^s to 
reduce Eqs. (|3|) and ( ^.ISj ) to 



1 



2e(x) 


K- 


/dv 




\dx 


d r 


Uv 


-0 


dx 


\dx 



= 0, 



= 0, 



(2.18) 
(2.19) 



with 9{x) e (— oOj+oo). The first self-consistency equa- 
tion Eq. (2. It) takes the form 



dxe2«^(-) = 8A, 



and 



Fc(A) — ng dx 



(d0/da 



2gKX. 



(2.20) 



(2.21) 



When deriving Eq. (2.21), we have taken into account 
Eq. (2.18), as well as the periodicity of 9'{x). 

Saddle-point equations Eqs. ( 2.18| ) and (2.1E) corre- 
spond to the global T-breaking perturbation as is set by 
Eqs. (2.1) and (pTsj). For the local perturbation, the lin- 
ear in (/J term disappears from the functional Eq. ( |2.5| ) 
leading to the replacement of the terms {dv/dx— 1)^ and 
dv /dx ~ 1 by th e terms (dv/dxY + 1 and dv/dx in Eqs. 
( 2.18 ) and ( 2.19| ) , r espec tively. It can be shown |jl^ that 
the modified Eq. (||l|) and the periodicity of v{x) re- 
sult in the condition dv/dx — which, in turn, leaves 
the only possibility to satisfy the modified Eq. ( ^.18 ) by 
a spatially homogeneous configuration of 9{x) at k = 1. 
This corresponds to the universal Wigner-Dyson regime 
(see also the discussion (i) a fter E q. (|2.27| )). 

Now we are back to Eqs. ( 2.18| ) - ( |2.2l| ) which can be 
brought to a more convenient form. Resolving Eq. (2.19), 



and making use of the periodicity of v{x) we obtain 

dxe-2»»=AA. (2.22) 



This is the second self-consistency equation. Further, in- 
troducing the new notations 



w{x) = ^» - e, 



(2.23a) 
(2.23b) 



In terms of solutions to the above saddle-point equations, 
the functional Fc{\) ~ F[Q\, 



^ ins 7 -^sp ; 0] 



iKspX reads 



dx 



) ( ^ ip] —Kip cosh 0{x) 



2gn\ 



(2.17) 



we rewrite the only remained saddle-point equation Eq. 
( ^l8| ) as 



d'^w 
dx"^ 



+ cr^ sinh ii;(.T) = 0, 



(2.24a) 
(2.24b) 



It should be supplemented by the two self-consistency 
equations which we combine to 



(2.25) 



Here, the periodicity of w'{x) has explicitly been taken 
into account. Also, 

F,(A) = ^ I dx ('^ + 2gn\. (2.26) 



Equation ( 2.24 ) is an equation of motion of a classical 
particle with a unit mass in the potential 



U{w) — CT^[coshi 



(2.27) 



w and X playing the role of coordinate and time, respec- 
tively. Due to a ring topology of our problem, we are 
interested in those solutions to Eq. ( p. 24 ) that describe 
the finite motion of the particle with the period T = 1. 
There are two kinds of solutions available. 

(i) The first solution is a trivial one, w{x) = 0. In this 
case, one concludes from Eq. ( 2.25| ) that t he sa ddle-point 
value oi Ksp is given by k = 1. Equation ( 2.26 ) then tells 
us that Fc{X) = F^^{X) ^ 2g\ for arbitrary A > 0. One 
can verify, by making use of Eq. ( 1.1b ) of Introduc- 
tion for /3 = 1, that this trivial solution appears to de- 
scribe (with the exponential accuracy) the distribution 
of the level curvatures in the universal, Wigner-Dyson 
regime. 

(ii) There is, however, the second solution which is 
spatially inhomogeneous. The important point here is 
that the finite motion in the potential U{w) cannot oc- 
cur with an arbitrary period. It is bounded from above, 
T{£) < Tinax = 27r/cr, with 



T{£) = \/2 



-{£) 



dw 



4K(fc) 



^£ - U{w) 



(l-fc2)i/4, 
(2.28) 
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Here, w± are the turning points of the finite motion, 
U{w±) = £, and £ is the integral of motion 



The parameter /c e (0, 1) is defined by 



K(A;) =F(7r/2,fc), with 



F(V^,fc)= / 



da 



sin^ a 



(2.29) 



(2.30) 



(2.31) 



37t 



2iz 



o 

LL 

'C3) 



g 



^/ 16g/(57t') 



-1 ,— inS/,^ V 

g (k) 




being the elfiptic integral of the first kind ||2(|| . 

The maximal period Tmax corresponds to the almost 
harmonic particle oscillations in a vicinity of the mini- 
mum w = of the potential U{w). Alternatively, for 
a given period T = \ (which we are interested in), the 
nontrivial period-1 solution exists only for a < Gc — 27r. 
It is given by 



1 



F(VXx),fc) 
2K(fc) '■ 



ipix) — arcsin ( — — \/l — fc^e'^f^) 



(2.32a) 
(2.32b) 



Equation ( 2.32 ) makes it possible to complete our com- 
putational program. To this end, we substitute Eq. 
( ^.32D into Eq. ( |2.25D to obtain 



E(fc) 



(2.33) 



Here, E(/s) = E(7r/2, k) is the elliptic integral of the sec- 
ond kind [Bel 



E(V',fc)= r daVl~k^ sin^a. (2.34) 

^0 



Further, we insert Eq. (2.32) into Eq. (2.26) and take 
into account Eqs. (|2!29|), (^.33|), (|^, ( |2.24bD as well 
as the identity 



cr = 4K(fc)(l-A:2)i/4 



(2.35) 



following from Eq. (2.28) taken at T = 1, to derive after 
some algebra the instanton solution Fc{X) = F™^{X) in 
the form 



F_;"^(A) = 47r5K(fc) [(2 - k^)K{k) - E(fc)] 



(2.36) 



FIG. 3. The function Fc{X) given by Eq. ( |2.39[ ). The 
dashed line depicts the Wigner-Dyson solution extended to 
the region A > Ac. A jump of the 'heat capacity' Cc(A) at the 
critical point A = Ac is shown in the inset. 

It remains to find the relationship between the initial 
parameter A of the Fourier space in the l.h.s. of Eq. 
( p. 36 ) and the auxiliary parameter k. This is done by 
combining Eqs. ( ^.24b| ), ( |2.25| ), ( ^.33] ) and ( |2.35D with 
the result 



2ti E3(fc) 
1 - P K(fc) ' 



(2.37) 



Let us analyze Eqs. (|2.36| ) and ( 2.37 ) which are man- 
ifestly diff'erent from Fj^^{\) = 25A with A > 0. It was 
already argued that the nontrivial solution Eq. ( 2.32 ) 
occu rs at a < ac — 2%. As is seen from Eqs. ( 2.35 ) 
and (2.37), this domain maps to < /c < 1 with kc — 0, 
and/or to A > A(0) = Ac, with 



A. 



TT 

Y 



(2.38) 



Exactly at the point A = Ac one observes that F^"''(Ac) = 
Fc^°(Ac), while at A > Ac the inequality ^^'""(A) < 
Fc^°(A) holds. [In particular, it follows from Eqs. ( |2.36| ) 
and ( p7| ) that |^ F^'^W ^ ngln \ as A > 1.] 
This means that for A > Ac it is the spatially inho- 
mogeneous, periodic solu tion Eq. ( 2.32 ) to the saddle- 
point equation Eq. ( 2.24 ) that minimizes the functional 
F[Q; K, ip] + iKX in Eq. (fj) as ^ 0. Therefore, we 
are forced to conclude that 



2gX for < A < Ac 
F^""(A) for A > Ac. 



(2.39) 



The result Eq. ( p. 39 ) supports the qualitative arguments 
given after Eq. (2.£). 
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A distinctive feature of Eq. ( 2.39| ), which holds as long 
as g ^ 1, is that the second derivative of FcW with re- 
spect to A exhibits a discontinuity exactly at the critical 
point A = Ac, in whose vicinity |A — Ad ^ Ac we derive 



^ = 2A-^e(A 

g 5 



(2.40) 



One can define the 'heat capacity' Cc(A) — —(PFc/dX^ 
which undergoes a jump at A = Ac: Cc(A) — for 
< A < Ac, and Cc(A) = IGg/dir^ for A = Ac + 0. 
This is the same discontinuity which is characteristic of 
the usual second-order phase transitions treated in the 
mean-field approximation, the factor g playing the role 
of the 'volume'. Equation ( 2.40 ) is an evidence to the 
second-order phase transition in the effective statistical 
mechanical model Eq. (I.IC) proposed in the Introduc- 
tion. 

The above discontinuity (revealed in the Fourier space 
A) has an important impact upon the distribution of the 
level curvatures. Being translated to the space of genuine 
curvature variable K, Eq. (2.40) implies appearance of 
the weak oscillations covering the main body of the dis- 
tribution function of the level curvatures, with a system- 



dependent amplitude of order 
period 



e ^ and the universal 



Ac TT^ 



(2.41) 

We stress that the value of 6K does not depend on mi- 



croscopic parameters of the Hamiltonian Eq. (2.1). Its 
value is only determined by the ring topology of the dis- 
ordered conductor as well as by the orthogonal symmetry 
of the problem weakly broken by the global external per- 
turbation in the form of a constant vector potential. 

We close this section by noticing that the above uni- 
versal oscillations do not appear for the local pertur- 
bation. In that case, neither the statistical mechanical 
model Eq. ( |1.10D can be introduced, nor are the non- 
trivial periodic configurations of the Q-matrix available 
in associated nonlinear tr-model p^ ]. This signals of a 
topological origin of the discovered oscillations. 



III. DISTRIBUTION OF LEVEL VELOCITIES AT 

P = 2 

In this section, we extend the previous analysis to 
study manifestations of the gauge invariance in the dis- 
tribution function of the single electron level velocities 
in quasi-onc-dimcnsional disordered ring where the time 
reversal symmetry has been broken from the very begin- 
ning so that the additionally applied small magnetic field 
respects the initial (unitary) symmetry of the problem. 
As is the case of the curvature distribution, the distribu- 
tion P^V) = A(^„ S{V - Vn)SiE - En{0))) of the level 
velocities 



A If 



(3.1) 



can be expressed through the parametric two-level cor- 
relation function R{lj, cj)) by means of the limiting proce- 
dure 



P^(y) = lim ifR {uj = VifiA, ip) , 

ip—>Q 



(3.2) 



as first proposed in Ref. pi] . 

We take the conventional route, and express the para- 
metric two-level correlation function in terms of the non- 
linear (T-model. The desired distribution of level veloci- 
ties is 

P,{V) = lim (^Re / VQp^Q] exp {-F[Q; V, ^]} . 

JqHx)=1s 

(3.3) 

Here, Pu[(5] is some preexponent which is irrelevant for 
our calculations performed with the exponential accu- 
racy. The free energy functional in Eq. (|3.3D reads 



(3.4) 



F[Q; V.^] = Y J dxStr {O^Q - i^[f, Q] 
iirVif 



dxStT{AQ). 



Despite the similarity of Eq. (3.4) to Eq. (2.5), the sym- 
metry properties JTot of the 8x8 supermatrix field Q are 
different in these two situations. [In particular, under as- 
sumption that the time reversal symmetry is completely 
destroyed, one has |||,|^ [Q,T3] = 0.] 

Evaluation of the Fourier transform of the distribution 
function of level velocities is the main goal of this section. 
In the instanton approximation used above one obtains 
(with the exponential accuracy) 

P^{X)^ I dV^P„(V^)e-^^^ocRee-(^['3"-^^=-"l+^^-^\ 



(3.5) 

where A is considered to be positive, A > 0, owing to the 
evenness of Pv{V). Let us proceed as follows. Due to 
the limit </? ^ to be implemented, we have to retain 
only noncompact degrees of freedom in the supermatrix 
Q which are given by its boson-boson sector. Simul- 
taneously, we omit the Grassmann entries, as discussed 
in the previous section. As a result, we replace 8x8 
supermatrix Q by conventional 4 x 4 mat rix in the 
parameterization of Eqs. ( 2.1C ) and ( 2.11 ), in which 



cos 6b — cosh 01 12, 
syhOb = i sinh0i]l2. 



(3.6a) 
(3.6b) 
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Here, </>, X € (0, 27r) and 9i G (0, +00) are spatially depen- 
dent, periodic functions of the variable x G (—1/2, +1/2). 
Substituting Qb{x) into Eq. (3^), and varying the func- 
tional F[Q\V,if\ -f iV\ with respect to the functions 
9i{x) — 0{x), 4>{x) and x{x)i we find the following set 
of the saddle-point equations that determine the instan- 
ton configuration Qiaa- 



smhO{x) 



dx' 



x{x) — ip I cosh6{x) 



= 0, 
(3.7) 



d 

dx 



dx' 



x{x) — f ) sinh^ 0(x) 



= 0. 



(3.8) 



Here, we have introduced k, = 2iT4p/g. Also, one has 
(j){x) = const. 

Variation over V yields the first self-consistency equa- 
tion 



TTtf I dx cos\v9{x) ~ A. 



(3.9) 



In terms of the above functions, the functional Fy{\) 
F[Q[ns; 14p, 0] + iVspX takes the form 



lim 



7rg 
2 



dx 



^] +sinh^0(.)f^ 
dx J \dx 



(f ) — K(p cosh0(a;) 



(3.10) 



To implement the limit cp ^ 0, we perform the following 
transformation, v{x) = x{x)/^ and 9(x) = 9{x) + ln(p, 
that results in two saddle-point equations 



dx'^ 4 



d_ 

dx 



dv ^ 
dx 



0, 



dv 
dx 



1 e 



(3.11) 



(3.12) 



with 9{x) e (—00,4-00). The first self-consistency equa- 
tion is 



TT / da;e^^(^) = 2A, 



while the functional F„(A) equals 



(3.13) 



(3.14) 



In deriving Eq. (3.14) we have used the periodicity of 
9'{x), and Eqs. ( plcf ), ^J)\ ) and ( |3l3| ). Second self- 
consistency equation is obtained from Eq. ( ^.12 ), 



(3.15) 



It is convenient to pass on to the new variables 
w{x) = ^9{x)+^, 



to reduce the saddle-point equation Eq. (3.11) to 



dx-' 

n 3 / K \ 2/3 



4 \Af 



(3.16a) 
(3.16b) 

(3.17a) 
(3.17b) 



The two self-consistency equations take the form 



1 



2A 



AA1/3k2/3 ^ 



AA2/3«;l/3. (3.18) 



Also, 



(3.19) 



Equation ( 3.17 ) describes a motion of a classical par- 
ticle of a unit mass in the potential 



2e2"'/3 + e"*^/^ - 3 



(3.20) 



w playing the part of the coordinate, and x being the 
time. Due to the Aharonov-Bohm geometry of the ini- 
tial problem, we are interested in the periodic solutions 
to Eq. dsTzl ), T = 1. 

(i) One such solution is a trivial one, w{x) — 0. In this 
case, we conclude from the two self-consistency equations 
Eq. ( p8| ) that Ari/3/«2/3 = 1 and 2\N^/^k^I^ = ^. It 
then follows that the saddle-point value of V^p is given by 
K — 2A/7r. Substituting the latter result into Eq. ( 3.19| ) 
we get F^{\) = F^^) = 9^lJ^^ aU A > 0. One 
can verify, by using [g8| Eq. ( 1.1a ) of Introduction, that 
the trivial solution w ~ corresponds to the universal 
Wigner-Dyson form of the distribution of the level ve- 
locities at /3 = 2. 

(ii) There is, however, a nontrivial period-1 solution 
to Eq. (3.17) provided A > Ac, Ac being some universal 
constant. Depending on the integral of motion 



2\dx 



3a' 



2e 



2w/3 



-4w/3 



(3.21) 
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the spatially inhomogeneous solution to Eq. (3.17) has 
the period 



J 1L 



dw 



_(£) ^£ - U{w) 
2VG 



c^\/7o - 7- 



K(fc). 



70 



7± 



•\/3 — tan(a/3) 
^/3 + tan(Q;/3) 
1 

2sin2/3(Q,/2) 
1 



2sin2/3(Q,/2) 



cos a = 1 — 2 1 



9^ 



1 + 2 cos 
1 — 2 cos 

3 



a ± TT 



(3.22) 



Here, w±(£) are t he turning points, U{w±) = £. Nota- 
tions in Eq. ( 3.22 ) are as follows: 



(3.23a) 
(3.23b) 
(3.23c) 
(3.23d) 



If we fix the period of a finite motion to be T = 1, the 
period"l solution (that exists only for a < (Tc = 27r) is 
given by 



1 _ F(^(x),fc) 

2 2K(fc) ' 



1 /e2«'(^)/3 - 



(3.24a) 
(3.24b) 



At this stage, we are able to find a closed analytical 
expression for Fy{X) defined by Eq. (3.19). To this end, 
we substitute the solution Eq. (3.24) into Eq. (3.18) to 
derive 



J dxe2™(^)/3 = 7+ + (7+-7_) 



E(fc) 



(1 - P)K(fc) 



- 1 



(3.25) 



Further, we insert Eq. (^.24|) in to Eq. ( [3.16|), a nd take 
into account Eqs. (|3^), (^.18|), (^Isj), ( |3.17b| ) as well 
as the identity 



2x/6 



V7o - 7- 



K(fc) 



(3.26) 



following from Eq. ( |3.22 ) taken at T = 1 , to obtain after 
some algebra the instanton solution Fy{\) = F^"'^(A) in 
the form 



Prix) 



70 - 7- 



-K^(fc) 



7- +70 - (7+ -7-) 



E(fc) 



(1 - fc2)K(fc) 



- 1 



Here, the parameter A in the l.h.s. of Eq. (3.27) is connected to the auxiliary parameter k € (0, 1) as 



2ttV2 
\/7o - 7- 



K(fc) 



7+ + (7+ - 7-) 



The ab o ve re lat ionsh ip is obtained by combining Eqs. 
( |3.17b| ), ( ^Isl) , ( |3^ ) and ( ^ ). 

Let us analyze the solution derived. It obviously dif- 
fers from F^^ = gAV27r with A > (that is related to 
a trivial case w{x) ~ 0). First, the instanton solution 
Eq. (3.27) arises at a < ac ~ 2tt that corresponds to 
A > A(0) = A„ 

27r2 
V3" 



A. 



(3.29) 



Ac, one observes that 



= F^°(Ac), while at A > Ac the inequality 
< F^^{X) is fulfilled. [In particular, one has 



Second, exactly at the point A 

11 F„(A) = 27r5ln^A as A > 1.] This means that 
for A > A c the spatially inhomogeneous, periodic solu- 
tion Eq. ( 3.24 ) to the saddle-point equation Eq. ( 3.17 ) 
leads to such configurations Qins which become energeti- 
cally more favourable as compared to the trivial solution 
w = 0. As a result, we conclude that 



i^.(A) = 



.gA2/27r 



for < A < Ac 
for A > Ac. 



(3.30) 



E(fc) 



(1 - fc2)K(fc) 



- 1 



(3.27) 



(3.28) 



The most important feature of Eq. (3.30), which is valid 
for g ^ 1, is that the second derivative of Fy{X) with 
respect to A exhibits a discontinuity exactly at the crit- 
ical point A = Ac, in whose vicinity |A — Ad ^ Ac the 
following expansion holds 



FyiX) 



A2 
2tt 
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e(A-Ac) 



'4 



(3.31) 



The 'heat capacity' Cv{X) = —(PF^/dX^ is seen to posses 
a jump at A = Ac: C„(Ac -t- 0) - C„(Ac - 0) = 18g/ll7r. 

The above discontinuity, similar to the one observed 
in usual second-order phase transitions, signals that the 
main body of the distribution function of the level veloc- 
ities gets dressed by weak oscillations with the amplitude 
of order . Their period is universal, 



Ac TT 



(3.32) 
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for it does not depend on microscopic parameters of the 
disordered conductor. The value SV is only determined 
by the unitary symmetry of the system and by the global 
type of the external perturbation applied. 



IV. RING TOPOLOGY: 
QUASI-ONE-DIMENSIONAL VERSUS 
WIGNER-DYSON LIMIT 



The analysis presented indicates that the universally 
oscillating corrections are of a topological origin: A ring 
geometry of a conductor subject to a constant vector po- 
tential (that represents a global perturbation and sets a 
global gauge invariance on the level of the microscopic 
description) is a necessary condition for their existence. 
However, the very fact of appearance of oscillations with 
a universal period, raises the question about actual im- 
portance of quasi-one-dimensionality for the above effect. 

Indeed, if we suppose that the ring topology (literally 
taken in its geometric sense, Fig. 1) is the only condi- 
tion required to observe the topological universality, it 
could be a plausible assumption that, as long as g ^ 1, 
the same oscillating features must persist in the distri- 
bution functions of level curvatures and level velocities 
in the ergodic regime. In the language of the random 
matrix theory, this is equivalent to the assumption that 
invariant random matrices closed in a ring (see below) 
should lead to the same results, provided the dimension- 
less conductance g is replaced by the effective conduc- 
tance g{N) linearly scaling with the matrix dimension N, 
g{N) ^ N. Below we argue that this is by no means the 
case. It turns out that quasi-one-dimensionality is an es- 
sential ingredient of the theory, being responsible for the 
forming the periodic boundary conditions in associated 
nonlinear a-model. It is these boundary conditions which 
ultimately led to identifying the universal critical point 
Ac in the Fourier transforms of the distribution functions 
in question and motivated our claim about existence of 
the second-order ph ase t ransition in the statistical me- 
chanical model Eq. ( LlOj ). 

To support this point, it is useful to develop a dis- 
crete (T-model approach to the periodic banded matrices 
defined by Eq. (L4). The advantage of this class of 
random matrices lies in the fact that, in the thermody- 
namic limit, they offer a possibility of probing the two 
physically different regimes depending on the behavior 
of the scaling parameter rj = b'^/N as TV — > oo. Namely, 
(i) If the band width b scales with the matrix dimen- 
sion iV as 6 = ^y^iVi/^, g^ 0{N°), the banded ran- 
dom matrices correspond to the quasi-one-dimensional 
physics {rj —^ g,); (ii) If the band width b scales lin- 
early with N, b ~ aN, < a < 1, the banded random 
matrices flow toward the ergodic, Wigner-Dyson regime 
(ji cx a^N — > oo). Below, only the main points of the 



derivation of nonlinear cr-model are sketched. For sim- 
plicity, the unitary symmetry is assumed; also, we set 

Let us represent the ensemble of periodic random ma- 



trices Eq. (1.4) in the form 



(4.1) 



where H^^^ and T?*-^-* are statistically independent, Her- 
mitian, Gaussian distributed matrices with the variances 



/^(l,2)^(1.2), ^ c j(l,2) 



(4.2) 



Here, J^u) describes the correlations between the matrix 
elements belonging to the band of the width 2b along the 
main diagonal, and J^n, is due to the correlations between 
the elements of the right upper and the left lower corners 
in Eq. (1^) (see Fig. 2). To make further consideration 
as clear as possible, we specify these correlation matrices 



J(2; =2e-^/''cosh[(Ai-i.)/6]. 



(4.3a) 
(4.3b) 



The choice Eq. ( p~^ ) enables us to treat the model of 
periodic banded matrices exactly. [In the case of con- 
ventional banded matrices, where only J^J appears, this 
fact has been used in Ref. pOt.] 

To derive a nonlinear cr-model, we follow the standard 
route |]l^,|^. The generating functional Z[J) reads 

Z[J)^ j dl'H] ^e-'*(^''*+'f^'")*^ g»*(£;+j+A(c^+«)/2)*^ 

(4.4) 

Here, is the source term, 5* is AN component super- 
vector, with = L = diag(+l, -1, +1, +1)r_a 
and A — diag(-f 1, H-1, — 1, — 1)r_a- Performing the av- 
eraging (...) over the ensembles of random matrices i?*^^^ 
and 77^^' , one obtains 



z{j) ^ / dm 



(4.5) 

where J^i, = J^],^ -I- J^l? , and = ^'^ (g) respects the 
symmetry = Lqj^L. The next s tep consists of decou- 
pling the interaction term in Eq. (4.5) by means of the 
Hubbard-Stratonovich transformation, followed by inte- 
grating out the supervector ^E*. Omitting the details, we 
show the result: 



N 

Z{J)= [ [^Pa^e^^EZ^.J^-^l-S'-^C^M--) 



N 



Y[Sdet-^E + J+^^A~aX (4.6) 
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The following observation is important here: Since the 
interaction term in the exponent of Eq. (4.5) depends 
only on that part of supermatrices which can be ex- 
panded in the eigenvectors 0^{fi) of the correlation ma- 
trix J^i, , the auxiliary supermatrices tr^ appearing in the 
course of the Hubbard-Stratonovich transformation and 



entering Eq. (4.6) have necessarily to belong to the same 
class of supermatrices specified by the expansion 



Co-Oo-(^), 



(4.7) 



Co- being the supermatrix coefficients. Equation (4.7) im- 
plies, in turn, that the supermatrices tr^ are not allowed 
to obey arbitrary boundary conditions. The boundary 
conditions are uniquely fixed by the structure of the cor- 
relation matrix J^,y (or, equivalently, by the structure of 



the inverse matrix [J ^]^,/). 

Although the generating functional in the form of Eq. 



(4.6) is only an intermediate (though exact) result in de- 
riving a nonlinear a-model, the conclusion made about 
the boundary conditions in a discrete cr-model is quite 
general: The constrained supermatrices Q^, = I4, 
that appear on the later stage of the saddle-point eval- 
uation of the iV-fold supermatrix integral in Eq. ([4.6|), 
must also belong to the same class, Eq. (4/7). Hence, 
the boundary conditions for constrained supermatrices 
are dictated by the structure of [J~ 



as weh. 



From now on, we are interested in the inverse matrix 
[J~^]tj.u- It possesses a remarkably simple tridiagonal 
structure supplemented by (generically) nonvanishing en- 
tries (1,7V) and (TV, 1): 



[J ]fJ.L' 



1 



1 - e-2/b 



(1 + e-^/'')2 - 
(1 + e-^/b f - e-2/b 
1 + e-^/" 



(l + e-A'/b)2 _e-2/b 



(4.8) 



At finite N, the entries (l,iV) and (A^, 1) are the hall- 
mark of the ring geometry incorporated in the random 
matrix Eq. ( [T^ ) at = 0. 

Let us analyze the eigenvectors Oa{p) of the inverse 
correlation matrix Eq. (4.8) for two different thermody- 
namic limits (i) and (ii) introduced prior to Eq. (4.1). 



conditions imposed on in nonlinear cr-model. [On the 
formal leve l, this is due to the (1, A^) and (A^, 1) entries 
in Eq. (O) that remain finite as A^ — > cx) and couple the 
fields (Ti and (Tjv in the exponent of Eq. (4.6).] 

(ii) In the limit b = aN , that corresponds to approach- 
ing the ergodic, Wigner-Dyson regime, the normalized 



(i) In the quasi-one-dimensional limit, b = gj N^/'^^ 



the normalized eigenvectors of [J' 
matrix 



are those of the 



eigenvecto rs o f [J 
with Eq. (EjI)) 



are those of the matrix (compare 



2<5„ 



{&ij..u+i + ^u,fi+i) ~ {5p^,i5^,i + 5^j^,n&u,n)- (4:. 13) 



25,. 



(<^p,i/+i + <^i/,^+i) — (<5^,i<^i/,w + '5iy,i'5^,7v), (4.9) Notice that the parameter a has dropped from Eq. 



which turns out to be independent of . The eigenvec- 
tors are found to be 



(4.13), and the entries 8^^i5u,n and 6i,,i5^j,^n have dis- 



1 



where 



{2TT/N)a for A^ even, 
(27r/A^) (o- - 1/2) for A^ odd. 



(4.10) 



(4.11) 



appeared, leading to destroying the coupling between cti 
and (Tat in t he exponent of Eq. (4.(;). The eigenvectors 
of Eq. (|4T|) are 



o.(m) 



cos[(w/A^)(/i-l/2)], (4.14) 



and 1 < cr < A^. 

It is seen from Eqs. (^7|), (|4.10| ) and ( [4.1l| ) that the 
supermatrices obey the periodic boundary conditions 



(4.12) 



where < ct < A^ - 1. 

An immediate consequence of the solution Eq. ( |4.14| ) is 
that the periodic boundary conditions for (and hence 
for Qfj) are not dictated by the above set of the eigenvec- 
tors Oct {p) ■ Contrary to the quasi-one-dimensional situa- 
tion, we deduce from Eqs. (47) and ( 4.14 ) the boundary 
conditions of the closed disordered sample which (in the 
continuum limit) read: 



as do the supermatrices of the corresponding nonlin- 
ear cr-model. Thus, we conclude that in the thermody- 
namic limit N 00 with the fixed parameter bN~^^'^, 
the geometric periodicity built in the matrix Hamilto- 
nian Eq. (11.4D is equivalent to the periodic boundary 



d . 



d_ 



0. 



(4.15) 



fi=N -1/2 



Equation (4.15) requires the current across the effective 
boundary of the sample to vanish. [Notice that Eq. 



12 



(4.15) is a close analog of the so-called 'trapping plane' 
boundary conditions arising in the radiative transfer the- 
ory [0.] Hence, we infer that in the thermodynamic 
limit N —> oo with the parameter bN~^ being fixed the 
connected geometry of the matrix Hamiltonian Eq. ( |1.4| ) 
is generically not equivalent to the periodic boundary 
conditions in the associated (nonlinear) tr-model. 

Equations (4.12) and ( 1.15[ ) are the main outcome of 
this Section. They show that in the thermodynamic limit 
N oo the closed geometry of the matrix Hamilto- 
nian Eq. ( |l.4| ) does not necessarily induce the period- 
icity of the supermatrix fields in the associated nonlinear 
cr-model. The criterion to distinguish between the two 
types of the boundary conditions follows from Eq. ( |4.S| ) 
that enables us to estimate the relative magnitude of the 



coupling entries <5^,i5i/,Af {Sv,iSi^,n) which equals be' 



-N/b 

provided 1 b N. The fulfilment of the inequality 
^g-JV/b ^ I leads to the periodic boundary conditions; 
vice versa, the inequality be~^^'' ^ 1 correponds to the 
boundary conditions of the closed sample. It was ob- 
served, that while such a periodicity is respected in the 
quasi-one-dimensional limit (i), it breaks down in the er- 
godic (Wigner-Dyson) limit (ii) being naturally replaced 
by the boundary conditions of the isolated disordered 
sample. This latter circumstance rules out the very pos- 
sibility to set up the problem of parametric level statistics 
with a gauge-type perturbation in the ergodic regime. 
For this reason, we have to conclude that the quasi-one- 
dimensionality is of vital importance for existence of the 
universal topological oscillations in the distribution func- 
tions of the level velocities and the level curvatures. 



V. CONCLUSIONS 



The consideration presented in this paper has been 
concentrated on two particular measures of parametric 
level statistics in quasi- one- dimensional conductors with 
a ring topology in the presence of a constant vector po- 
tential. They are the distribution functions of the level 
curvatures (/3 = 1) and the level velocities (/3 = 2 and 
4). We have used the instanton formulation of the re- 
duced nonlinear cr-model to evaluate nonperturbative, in 
inverse Thouless conductance g~^, corrections to the uni- 
versal random-matrix-theory distribution functions Eqs. 
( ITTal) and ( p^ . Our main findings are rather startling, 
especially in the light of the previous studies on the sub- 
ject, which have created a common belief that small but 
finite value of inevitably results in nonuniversal cor- 
rections beyond RMT. Contrary to these expectations, 
we have demonstrated that in quasi-one-dimensional dis- 
ordered conductors with g ^ 1 there exists another type 
of corrections which, to a large extent, are universal. The 
corrections are of topological origin and manifest them- 
selves in the form of the weak oscillations on the main 



bodies of the distribution functions of the level curva- 
tures and the level velocities. While the magnitude of 
the topological oscillations is system dependent (~ e^^), 
their period appears to be independent of microscopic pa- 
rameters of conductor, and hence universal. It is entirely 
determined by the global symmetries of the Hamiltonian 
before and after the perturbation was applied. We have 
predicted the period of the oscillations to be 4/7r^ for 
the distribution of level curvatures at /? 1 , and VS/tt 
for the distribution of level velocities at /3 = 2 and 4. It 
should be stressed that no rescaling is needed to establish 
the universality of the topological oscillations: Their pe- 
riod is universal and parameter independent in genuine 
curvature and velocity variables. 
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APPENDIX: DISTRIBUTION OF LEVEL 
VELOCITIES AT /3 = 4 

The appendix aims to demonstrate that the conclu- 
sions drawn for the distribution function of the level ve- 
locities dX (3 — 2 remain valid for /3 = 4. The starting 
point of our analysis is Eqs. (3.5) and (3.4), in which the 
parameterization of the supermatrix Q{x) must corre- 
spond to the symplectic symmetry pO| . Due to the limit 
(/? — > to be taken, we keep only boson-boson sector in 
Q; also, the exponential accuracy of our calculations al- 
lows one to eliminate the Grass mann entries there. We, 
therefore, arrive at Eq. ( 2.1(]| ), with 6b given by Eq. 
the matrix U{x) is to be replaced by 



U{x) 



I2 

F{x) 



with 



F{x) = 



D{x) 



(A.l) 



(A.2) 



X [{D{x) — 2)l2 + 2imc7z + 2im[aj. + 2im'(ay] , 
D{x) = 1 + m^{x) -1- mf{x) + mf{x). (A.3) 



Here, 



~i 

1 



(A.4) 



is the y-Pauli matrix [a^ and were specified 
by Eqs. ( p. 13 ) and (p.^)), and the functions 
'm{x)^mi{x),m'({x) G (—00, +00). In the above param- 
eterization, the free energy F„(A, ip) = F[Q; V, ip] -t- iVX 
is 



13 



cosh6'(x) + ip'^ sinh^ 6{x) > + 27rg / dx 



sinh 



dm 
dx 



/ dm'i 
V dx 



dx I \ dx dx dx dx 



(A.5) 



In order to simplify the formulas, we perform the limit — > in Eq. ( |A.5| ) prior to taking the variational derivatives. 
This is done by absorbing tp into 9{x) = 0{x) — hiip, followed by introducing the new functions h{x) = — 2m(x)/<y9, 
h!i(x) = —2m'i{x)/Lp and h'l{x) — —2m'({x)/ip. The limit ip ^ results in 



dh ^ 

dx 



dh[ 
dx 



dh'i 
dx 



iV{\-- I da;e^(^' ) . (A.6) 



Evaluation of the instanton action is by now standard. Taking the functional derivative of Eq. 
the functions entering Fy{X,0), we obtain the following set of the saddle-point equations: 



with respect to 



d'e 1 



dx^ 



dh ^ 
dx 



dh[ 

dx 



dh'l 
dx 



0, 



(A.7) 



d 

dx 



— -\\ e2S(^) 
dx 



d_ 

dx 
d_ 
dx 



dx 
dh'l 
dx 



= 0, 
= 0. 



(A.8) 
(A.9) 
(A.IO) 



Here, k = 2iVsp/g. Variation over V yields the (first) 
self-consistency equation 



TT / dxe^^'='> = 2A. 



(A.ll) 



Let us first consider Eq. (A.9). It is equivalent to 

'^^1 _ (J g-2e(x) 
dx 



(A.12) 



Integrating it over x G (—1/2, -1-1/2), and taking into ac- 
count the periodicity of /i'i(x) we conclude that C = 0, 
so that 



dh{ 
dx 



0. 



Analogously, we conclude from Eq. (A.IO) that 



dh'l 

dx 



0. 



(A.13) 



(A.14) 



Equations ( A.13| ) and ( A.14|) enable us to write the 
instanton free energy F^{\ ) — f [Qinsi ^4p,0] -I- iVs^X ob- 
tained by minimizing Eq. ([A.6|) in the form 




dx I 4 \ dx 



(A. 15) 



where the functions 9{x) and h{x) obey Eq. 
the equation 



dx'^ 4 



, dh 

K- ( — -1 

dx 



= 



(A.16) 



that follows from Eqs. ( A.7|) ( |A.13| ) andjATj. 

Applying now Eq. ( |a.16D to Eq.^ (|A.15D , and tak- 
ing into account the peri odicity of 0'{x) and the self- 
consistency condition Eq. ( A.ll ), we come down to 



'Kg 
~2 



dx 



dx I 



-gnX. 



(A.17) 



I) and At t his po int we notice that E qs. (A.17|), dA.16|), ( |A.8| ) 



and jA.llD coincide with Eqs. ( |3.14| ), ( p.llD , ( |3.12D and 
( |3.13|) , respectively. This proves that the universal, non- 
perturbative in g~^, corrections to the distribution of the 
level velocities at /3 = 4 coincide with those found with 
the exponential accuracy for (3 = 2. 
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